Mass dependence of nuclear shadowing at small Bjorken-x from diffractive scattering 
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I. INTRODUCTION 

Calculations of various cross sections in relativistic 
nucleus-nucleus collision physics in the framework of per- 
turbative QCD (pQCD) require nuclear parton distribu- 
tion functions (nPDFs). The nPDFs are usually pre- 
sented in a parameterized form, as a product of free- 
nucleon PDFs and a "shadowing function" (see e.g. pj), 
which depends on x and momentum transfer Q 2 . Exper- 
imentally extracted shadowing functions are limited by 
the kinematic range of the data, and fits to these data do 
not provide the shadowing function in the whole requisite 
(x,Q 2 ) region for all calculations. 

Another approach is to utilize Gribov theory[2j, which 
relates shadowing to diffraction. For the deuteron, the 
relationship involves the interaction of the diffractively 
scattered (virtual) photon with only two nucleons, but 
for heavy nuclei triple and higher-order scattering may 
be important and needs to be included in the formal- 
ism. This leads to some model dependence. For nuclear 
collision applications, nPDF uncertainties contribute sig- 
nificantly to the overall uncertainties of calculated quan- 
tities. Thus, an accurate determination of nPDFs is of 
paramount significance. 

It is well-established experimentally that for small val- 
ues of the Bjorken variable (x < 0.1), the nuclear struc- 
ture functions F^ are significantly reduced compared to 
A times the free nucleon structure function , where 
A denotes the mass number of a specific nucleus. Equiv- 
alently, the virtual photon-nucleus cross section is less 
than A times the one for free nucleons, a~ ( *A < Aa 7 *N. 
This phenomenon became known as nuclear shadowing 
(in a strict sense). The same behavior is observed for real 
photons at sufficiently high energies [y > 3 GeV). Shad- 
owing constitutes the most pronounced nuclear effect in 
lepton-nucleus deep inelastic scattering (DIS). (Note that 
in a generalized sense all modifications of the structure 
function in the entire range of x, including regions where 
F2 > A*F^ , are sometimes referred to as "shadowing" .) 

The shadowing ratio can be defined as F^/(A * FJ?). 
Since the virtual photon-nucleus cross section is propor- 
tional to i*2j the shadowing ratio can be expressed al- 
ternatively as o^-aI^A * a-y-N). In this work we employ 
Gribov theorv. including the real Dart of the diffractive 



scattering amplitude, to calculate the shadowing ratio at 
very small Bjorken-x. We compare to experimental re- 
sults from the NMC[1 H| and E665(l If experiments. 
These experimental data are all at small Q 2 ; we thus use 
the information from diffractive scattering of real pho- 
tons (Q 2 = GeV 2 ). 

The paper is organized as follows: in Sec. [H] we review 
the basic formalism of Gribov theory as applied to shad- 
owing in the small Bjorken-x regime. This is essentially 
the same treatment as described in Ref. Q , where the ba- 
sic object is the diffractive structure function. Here the 
basic quantity is the differential diffractive cross section. 
Sec. IIIII deals with diffractive production and the diffrac- 
tive dissociation cross section. We present the results of 
our calculation in Sec. IIVI and conclude in Sec. |Vj 



II. DIFFRACTION AND NUCLEAR 
SHADOWING 

A. Shadowing Ratio 

Several natural length scales may become important 
in scattering off of a nucleus of mass number A. Most 
important for the present discussion is the average inter- 
nucleon separation, d, of the order of 2 fm. A high energy 
(virtual) photon scattered from the system can experi- 
ence nuclear effects in two distinct ways|8|: 

i) incoherent scattering from A nucleons with mod- 
ifications to the structure functions due to many- 
body effects in the nuclear medium; 

ii) coherent scattering processes involving more than 
one nucleon at a time. 

The latter effects occur when hadronic excitations (or 
fluctuations) produced by the high-energy photon propa- 
gate over distances (in the laboratory frame) comparable 
to or larger than the characteristic length scale d ~ 2 fm. 
Shadowing can be understood in terms of coherent scat- 
tering on more than one nucleon. Incoherent scattering 
occurs primarily in the range 0.1 < x < 1, while strong 
coherence effects arc manifest at x < 0.1. 
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The (virtual) photon-nucleus cross section can be sep- 
arated into a part which accounts for the incoherent scat- 
tering from individual nucleons, and a correction (shad- 
owing correction) from the coherent interaction with sev- 
eral nucleons: 



Z <j~, 



(A - Z) a 



7*71 



(i) 



The single scattering part is the incoherent sum of 
photon-nucleon cross sections, where Z is the nuclear 
charge number, and er 7 * p and er 7 * n are the photon-proton 
and photon-neutron cross sections, respectively. The 
multiple scattering correction can be viewed as expanded 
in terms of the number of nucleons in the target involved 
in the coherent scattering (n > 2). The leading contri- 
bution to nuclear shadowing comes from double scatter- 
ing, since the probability that the propagating hadronic 
excitation coherently interacts with several nucleons de- 
creases with the number of nucleons. For the deuteron, 
only the double-scattering term [n = 2) is present. 

We define the shadowing ratio as 
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7*n 
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Thus, the evaluation of the shadowing correction, da^Ai 
is central to the calculation of the shadowing ratio. In the 
next section we utilize the Gribov theory in a generalized 
form to determine Sa 7 *A- 



Shadowing Correction From Generalized 
Gribov Theory 



The Gribov theory highlights the deep connection be- 
tween nuclear shadowing and diffraction. The original 
formulation relates the correction to the cross section in 
hadron-deuteron scattering to diffraction off of a nucleon. 
For the deuteron (and other light nuclei where double 
scattering is dominant), the shadowing correction can be 
evaluated without any ambiguity associated with the in- 
volvement of more than two target nucleons. For heavier 
nuclei, one needs to go beyond double scattering. Mul- 
tiple scattering with n > 2 nucleons may contribute sig- 
nificantly to the shadowing correction (depending on the 
kinematic regime). These contributions are not model- 
independent; therefore the evaluation of the shadowing 
correction using Gribov formalism depends to some ex- 
tent on the model for the n > 2 multiple scattering. 

The original formulation by Gribov neglects the real 
part of the diffractive scattering amplitude. Here we take 
this into account and denote by 77 the ratio of the real to 
imaginary parts of the diffractive scattering amplitude. 
In the generalized form incorporating the real part, the 
shadowing correction at the level of double scattering is 




FIG. 1: Double scattering contribution to nuclear DIS. 



given by 
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Here a^i s N is the photon-nucleon diffractive cross sec- 
tion and z represents the beam direction. The coher- 
ence length, A, is given by A = 2v/M x for real pho- 
tons. As illustrated in Fig. ^ a diffractive state with 
invariant mass Mx is produced in the interaction of the 
photon with a nucleon located at position (b, Z\) in the 
target. The hadronic excitation is assumed to propa- 
gate at fixed impact parameter b and to interact with a 
second nucleon at z%. The probability to find two nu- 
cleons in the target at the same impact parameter is 
described by the two-body density p A (b,zi; 6,^2) nor- 
malized as J d 3 r d 3 r' p A (r,f) — A 2 . The phase factor, 
exp{i[(zi — Z2)/A]} in eq. @ implies that only diffrac- 
tively excited hadrons with a longitudinal propagation 
length larger than the average nucleon-nucleon distance 
in the target, A > d ~ 2fm, can contribute significantly 
to double scattering. The limits of integration define the 
kinematically permitted range of diffractive excitations, 
with their invariant mass Mx above the two-pion produc- 
tion threshold and limited by the center-of-mass energy 
W = y/s of the scattering process. 

In order to approximate the two-body density 
Pa (b, Zi',b,Z2), we note that short-range nucleon-nucleon 
correlations are relevant in nuclei when zi — z\ is compa- 
rable to the range of the short-range repulsive part of the 
nucleon-nucleon force, i.e. for distances ^ 0.4 fm. How- 
ever, shadowing is negligible in this case. Short-range 
correlations are therefore not important in the shadowing 

domain, and the target can be considered as an ensemble 

(2) —i -* -* 

of independent nucleons with p A (r,r ) ~ PAv)pAv )j 

where pa is the nuclear one-body density. 

With increasing photon energies or decreasing x down 
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to x <C 0.1, the longitudinal propagation length of diffrac- 
tively excited hadrons rises and eventually reaches nu- 
clear dimensions. Thus, for heavy nuclei interactions of 
the excited hadronic state with several nucleons in the 
target become important and should be accounted for. 
Following we introduce an attenuation factor with an 
effective hadron-nucleon cross section, a c g . The shadow- 
ing correction can thus be written as 
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The effective hadron-nucleon cross section, a c s in eq. (@J 
is defined as 
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where er 7 jv is the photon-nucleon cross section. The de- 
tails of this approach and the approximations inherent in 
the definition of a e s are treated thoroughly in Q. For 
vector mesons as the intermediate hadronic excitations, 
we take ctvn as er c ff in the attenuation factor in cq. 0}, 
where uvn is the vector meson- nucleon scattering cross 
section. 

When x <C 0.1, it is a good approximation to ignore 
the phase factor, expi[(zi — z 2 )/A], in eq. Q, and (using 
Leibnitz's rule) the integrals over z\ and z 2 can be carried 
out explicitly. This leads to a simplified form, 

2(1 - l/A)a lN 



5 a 



~f*A 



Coff 




exp{-LT(b)} - l + LT(b) 



(6) 



where L = (A/2) (l-ir))a cS andT(fe) = JTr^dz p A (b, z), 
the usual Glauber thickness function. However, we use 
the full expression, eq. 10} , in our calculation. 

Eq. Q gives the shadowing correction in terms of 
CTeff- This effective cross section involves the differen- 
tial diffractive dissociation cross section at small t, see 
eq. (0 . The treatment of this diffractive cross section is 
the subject of the next section. 



III. DIFFRACTIVE DISSOCIATION 
A. Diffractive production 




FIG. 2: Diffractive scattering from a proton. 



dissociate, and it remains intact during the process. The 
photon, on the other hand, dissociates into a hadronic 
final state X, which is well separated in rapidity from 
the proton, 
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Such diffractive processes are important at small momen- 
tum transfer, with cross sections which decrease expo- 
nentially with the squared four-momentum transfer. In 
general they exhibit a weak energy dependence. 
Diffractive dissociation of real photons, 



~f + N ^ X + N , 



(8) 



has been studied in both fixed target and collider exper- 
iments. Experiments were carried out at Fermi National 
Laboratory (FNAL) at average photon-proton center of 
mass energies of W ~ 12.9 GeV and W ~ 15.3 GeV0. 
Diffractive states with an invariant mass squared of up 
to M\ ~ 18 GeV 2 were produced. This experiment 
measured the diffractive dissociation cross section differ- 
ential in both, the invariant mass Mj and the squared 
four-momentum transfer t. Experiments at the Hadron- 
Electron Ring Accelerator (HERA)[H [H [H [3 [l| 
were carried out at average energies W ~ 187 GeV and 
W ~ 231 GeV. Diffractive states with mass M x < 30 
GeV were produced. Unlike the FNAL experiment, only 
d<j^l s N /dMx was measured due to poor resolution in t. 

As mentioned earlier, the available experimental data 
on shadowing at small x (x ~ 10~ 4 ) are all at small 
Q 2 (Q 2 < 1 GeV 2 ). At such small virtualities the pho- 
tons can be considered quasi-real, and it is thus not a 
bad approximation to regard them as real photons with 
Q 2 = GeV 2 . The center-of-mass energies are also low: 
W ^ 15 GeV for the NMC and W ~ 25 GeV for the 
E665 measurements. These energies are comparable to 
the photon-proton center of mass energies at FNAL. For 
these reasons one can use the information from diffractive 
scattering of real photons at FNAL to calculate the shad- 
owing ratio in the kinematic range accessible at NMC and 
E665. 



B. Diffractive dissociation cross section 



Consider the single diffractive scattering of a (virtual) 
photon off of a proton (see Fig. 15}. The proton does not 



It is natural to divide the diffractive dissociation cross 
section data as a function of M\ into the region of the 
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with H v (Mjf ) the vector meson part of the photon spec- 
tral function which is given by 
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M v 2 [GeV 2 ] 



FIG. 3: (Color Online) Triple-pomeron fit to the FNAL 
dataHl , including the p' resonances in the continuum, as in 
Ref. [Li. 



low-mass vector mesons (p, u>, and <fi) and a high-mass 
continuum, with a matching point of M\ ~ 2.25 GeV 2 , 
as in Ref. |l5j . However, in this case a triple-pomeron 
fit to the continuum data (as described in Sec. IIII B 3|) 
misses the p' resonances as can be seen in Fig. 01 

The analysis by the HI collaboration^^ divides the 
HERA photoproduction data into effectively three inter- 
vals in M\. We follow this strategy in the present paper, 
taking the first interval (0.16 — 1.58) GeV 2 to contain the 
region of the low-mass vector mesons (p, to and <ft). The 
second interval (1.58 — 4.0 GeV 2 ) covers the p' resonance 
region. The third interval (M% > 4.0 GeV 2 ) is that 
of the high-mass continuum. The differential diffractive 
cross section is thus written as a sum over contributions 
from these three mass intervals, 
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(9) 
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In the following, we briefly summarize the various ap- 
proximations applied in the three regions. 



1. Low-mass vector mesons 

Generalized vector meson dominance (VMD) 0] can 
be utilized to describe the contribution of the low-mass 
vector mesons to the differential diffractive cross section, 
i.e. the first term on the right-hand side of eq. ©: 
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In eq. (|10p. cry at is the vector meson- nucleon cross sec- 
tion and e 2 /47r = 1/137 is the fine structure constant. 
The us and cf> mesons are narrow and thus well approxi- 
mated by delta functions. One can write their contribu- 
tion to the photon spectral function as 



n v (M 2 ) 



^ ) S (M 
9v 



V = W ,0,(12) 



where my and <?y, (V — u>,(f>) are the mass and the 
coupling constant of the u> and <fr mesons, respectively. 

The p-meson, unlike the u> and <f> mesons, has a large 
width due to its strong couplin g^ to two-pion states. We 
have followed the approach in |l5| and taken this into 
account through the 7r + 7r _ part of the photon spectral 
function: 
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(13) 



where is the mass of the pion and JWj = M nv is 
the invariant mass of the 7r + 7r _ pair. The pion form 
factor, Ftt is taken from Ref. A full discussion is 

given in [l5|]. We compared the result from the delta 
function approximation to this more exact calculation 
and found that taking into account the width of the p- 
meson increases the differential diffractive cross section 
by 10%. 

The vector meson-nucleon cross section in eq. <|10[) has 
an energy dependence of the form 

a V N -W 2 ^^- 1 ) =W 2t (14) 

where a w (t = 0) = 1+e is the soft pomeron intercept [l8j. 

2. Region of the p' resonances 

The p 1 resonance region contains the p(1450) and 
p(1700) mesons. These resonances were formerly clas- 
sified as the p(1600) I 11 R e f- [13 the p' resonance 
region is lumped together with the high-mass continuum. 
As can be seen in Fig. |3 the FNAL data show an en- 
hancement in this region. We treat this enhancement in 
terms of an average p' resonance, correspondingto the 
earlier classification of p(1600), as done in Ref. 0. We 
use the available information on the p(1600) from Ref. 
(l(j| in a VMD-type calculation to evaluate the contribu- 
tion from this region. The average p' resonance should 
have a finite width, but encouraged by the fact that a 



5 



delta function in the case of the p gives a good approxi- 
mation to the full-width result, we employ a narrow res- 
onance approximation for the p(1600). Thus, for the sec- 
ond term of © we have 



dM x dt 



tfsO 
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and V = ,9(1600). 
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3. High-mass continuum 

A full treatment of both the FNAL data and the HERA 
data in this region has been carried out by the HI Col- 
laboration in Ref. [10j], using the triple- Regge model of 
photon dissociation. For simplicity we consider only the 
triple-pomeron term, ignoring the sub-leading reggeons 
and interference between the pomeron and reggeons. In 
this case the differential dissociation cross section can be 
written as 



d 2 a 
AM 2 At 
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where so is the hadronic mass scale and B(W 2 ,M^ C ) = 
2b p jp + + 2a' w lnW 2 /M 2 . Here b p jp is the proton- 

pomeron slope parameter, 6p]pip the triple pomeron slope 
parameter, and a' w is the slope of the pomeron trajectory. 
We use the values in [ljj f° r these parameters, with s = 
1 GeV 2 . 

In order to have a consistent treatment of the different 
contributions, we take 0^(0) = 1.068 ±0.0492 from [To] . 
This agrees within error with the soft pomeron intercept 
in Ref. [l| (a w (0) ~ 1.081). Thus the only free pa- 
rameter is the triple-pomeron coupling Gppp(0). We fit 
the FNAL data at the average energy W — 14.3 GeV and 
t = -0.05 GeV to cq. The fit gives Gipipip(O) = 11.4 
Atb/GeV 2 ( 1+£ ). The quality of our fit is shown in Fig. 0] 



IV. RESULTS 

The treatment outlined in the last two sections is now 
applied to calculate the shadowing correction, and hence 
the shadowing ratio. The basic equation is eq. JBJ, which 
involves the ratio of the real to imaginary amplitudes ry, 
the photon- nucleon cross section cr 7 jv, the nuclear density 
Pa, and the effective cross section a c s- 

We use the e nerg y- independent 77 's for the vector 
mesons from Ref. IIm. For both p and ui mesons, ri takes 



M x [GeV'l 



FIG. 4: (Color Online) Triple pomeron fit to the FNAL0 
data in the nonresonant continuum (Mx > 4 GeV 2 ). 



values between and —0.3. Here we take r\ p = r\ u = —0.2 
in accordance with Ref. [l6^ . The results of our calcula- 
tion are not very sensitive to the precise values of rj p ^ . 
For the <p meson, we take -q^ = 0.13 [2£j. For lack of 
information, we take ^ p (i6nn) = 0. For the high-mass 
continuum, we follow Ref. [7j and define rj con t as 



Vcont = ^ ("ip (°) - !) » 



(18) 



using the result of Gribov and Migdal [2l|. 

The small difference between the photon-proton cross 
section a JP and the photon-neutron cross section a 
is neglected in this study. We 
Landshoff parameterization of cr 7P 
photon-nucleon cross section ct 7 jv F° r the nuclear den- 
sities three-parameter Fermi (3pF) distributions are ap- 
plied: 
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the Donnachie- 
as the generic 
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(19) 



with the parameter values taken from Ref. |22| . For mass 
numbers A < 20 a harmonic oscillator (HO) density dis- 
tribution may be more appropriate than the 3pF dis- 
tribution. For uniformity, we use the 3pF distributions 
for the whole mass range in light of the fact that un- 
certainties associated with other parameters are at least 
comparable. 

We carried out calculations at W = 15 and 25 GeV, 
the approximate energies of the NMC and E655 experi- 
ments, respectively. The results are displayed in Fig. 
together with the experimental data points. At very 
small x (x ~ 10 -4 ), NMC has two data points, corre- 
sponding to 6 Li and 12 C. The E665 experiment has four 
data points: 12 C, 40 Ca, 131 Xe, and 208 Pb. In view of 
the large error bars of the data, the calculation can be 
said to describe the experimental information in the en- 
tire mass range at E665 energies, and also for the two 
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FIG. 5: (Color Online) Shadowing ratio calculated at W — 
15 GeV (crosses) and W = 25 GeV (stars). Data are from 
the NMC (shaded circles) and E665 (shaded boxes) collabo- 
rations. The NMC point corresponding to 12 C is displaced 
slightly for better visibility. The dashed lines are energy- 
dependent fits according to 12011 as described in the text. 



A as W increases. This seems to suggest that the mass 
dependence of the calculated shadowing ratio at large 
A and increasing W is not as simple as in eq. 
A five-parameter fourth-degree polynomial with energy- 
dependent coefficients gives a good fit for the entire mass 
range and at all energies considered. However, we pre- 
fer the simple physical picture of the two-parameter fit, 
which is adequate considering the experimental and the- 
oretical error bars. 

The uncertainties of our calculation are mostly related 
to the various parameterizations of the diffractive dissoci- 
ation cross section. The delta function parameterization 
for the w and <f> mesons should be satisfactory, and the 
width of the p meson has been taken into account. Re- 
finements of the spectral function (|13[1 are possible, and 
improvements of the treatment of the p' resonance re- 
gion is also left for future work. The uncertainties in the 
continuum are associated with the neglect of sub-leading 
reggeons and of interference terms. Furthermore, the 
use of an effective scattering cross section to account for 
multiple scattering is an approximation, as is using real- 
photon information {Q 2 — 0) at small, but non- vanishing 
Q 2 . We estimate the overall uncertainty of our calculated 
results to be in the 20% range. 



available data points at NMC energies. For small A the 
shadowing ratio decreases rapidly with A, while for large 
A the decrease is more gradual. The calculated result 
also falls steeply with increasing W. 

The dashed lines in Fig. [S] represent a two-parameter 
fit of the standard form 

n s A = foA^- 1 (20) 

to the calculated results, with energy-dependent param- 
eters /?o and 0\ . In order to determine the energy depen- 
dence of the fit parameters we calculated the shadowing 
ratio for W in the range 10 < W < 30 GeV. The energy 
dependence of the fit coefficients can be described as 

0o = 0.720 + 0.1 18 \n(W) (21) 

and 

0! = 1.143- 0.075 bx(W) , (22) 

where W is the center-of-mass energy in GeV. 

The fit does very well for the entire mass range for 
low W and deviates from the calculated result at large 



V. CONCLUSION 

We have calculated the shadowing ratio at very small 
Bjorken-x for nuclei in the range 3 < A < 239 using Gri- 
bov theory. We included the effect of the real part of 
the diffractive scattering amplitude. The photon diffrac- 
tive dissociation cross section, which serves as an input 
to our calculation, was parameterized as a function of 
the invariant mass of the diffractively produced hadronic 
excitation using vector meson dominance and Regge the- 
ory in three mass intervals: low-mass vector mesons, p' 
resonances, and continuum. The parameters needed are 
taken from earlier studies that fit experimental data. 

It is found that the calculated shadowing ratio de- 
creases with mass number first rapidly and then more 
slowly. To be able to compare to NMC/E665 data, the 
calculations are all at low center-of-mass energy, and the 
decrease with A is stronger as the center-of-mass energy 
increases. We find that Gribov theory gives a reasonable 
estimate of the mass dependence of nuclear shadowing at 
small Bjorken-a;. 
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